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Abstract. We show that electronic wave functions ip of atoms and molecules 
have a representation ?/) = J^<j), where is an explicit universal factor, locally 
Lipschitz, and independent of the eigenvalue and the solution ip itself, and 
has locally bounded second derivatives. This representation turns out to 
be optimal as can already be demonstrated with the help of hydrogenic wave 
functions. The proofs of these results are, in an essential way, based on a new 
elliptic regularity result which is of independent interest. Some identities that 
can be interpreted as cusp conditions for second order derivatives of i/i are 
derived. 



1. Introduction 

1.1. Motivation and results. The non-relativistic quantum mechanical Hamil- 
tonian of an A^-electron molecule with L fixed nuclei is given by 

HnA^, Z) = -a + V{X, Z) + f/(X, Z), 

where V, the Coulombic potential, is given by 

N L 



(1.1) V^ViX,Z) = -}^}_^ 



,Xi. ~ Xi I ^-^ \xj 

j = l k=l ' ^ ' l<i<j<N ' ' 



and the internuclear repulsion U by 



l<k<e<L ^'^'^ '^^^ 

The latter is merely an additive term that will be neglected in the sequel and we 
will henceforth consider 

(1.2) H = Hn,l{X,Z)-U{X,Z). 

Above, X — {xi, X2, ■ ■ ■ , xn) G K'^^ denotes the positions of the TV electrons, with 
Xj = {xj,i,Xj^2,Xj^3) E the position of the j*'^ electron. The positions of the 
L nuclei with the postive charges Z = (Zi, Z2 ■ ■ ■ , Zl) g R;^ are denoted by X = 
{Xi,X2, Xl) e R^^ where Xu = (^/c,i, ^fe,2, ^fe.a) € is the (fixed) position 
of the A;*'' nucleus with charge Z^, and it is assumed that Xf^ ^ Xk for (. ^ k. 
The Laplacian corresponding to the j*^ electron is Aj = X]j=i ax ^^"^ 
Laplacian on R'^^ is given by A = X^jLi ^j- We also introduce the 3A^-dimensional 
gradient by V = (Vi, . . . , Vjv)- 
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The operator H is selfadjoint on L'^{M.^^) with operator domain V{H) = W'^^'^{M.^^) 
|14j . and it depends parametricaUy on X and Z. In the case of an A^-electron atom 
with (one) nucleus of charge Z fixed at the origin G K'^, (|1.2|) becomes 



Generations of chemists and physicists have devoted a good part of their research 
to the analysis of various problems related to ffAr,L(X, Z). Most of the present 
day understanding of atoms and molecules is based on the analysis of problems 
directly related to this operator, see any textbook in atomic and molecular quantum 
mechanics. 

One of the central problems is the eigenvalue problem 



Since the electrons are Fermions the A'^-electron wave function ip has to satisfy the 
Pauli Principle. This can be achieved in a spinless formulation by requiring that ip 
transforms according to certain irreducible representations of the symmetric group 
6^ . Our present work will not require any symmetry assumptions on ip. More 
precisely, we will consider local properties of distributional solutions (locally L^) in 
a domain C M.^^ to Hip — Eip where E can be any real number. 

Within mathematics and mathematical physics Schrodinger operators as (|1.2|) 
are studied mostly from an operator theoretical point of view, see the textbooks 
PP, [HI, [IT, and as weU as the recent survey pOl . 

The PDE-aspects of (|1.4f) have been studied in relatively few works. We first 
note the following: Let S(X) denote the set of points in R'^^ where the potential 
V defined in is singular. The function V is real analytic in 'E?^ \ S(X) and 
hence by classical results (see ^2 Section 7.5, pp. 177-180]), so is ip. 

Therefore a basic question is how to characterize the effect of the singularities 
of V on the local behaviour of a solution ip of (|1.4|) . 

In 1957 Kato showed that a solution ip satisfying H1.4|l is continuous in all 
of M?'^ with locally bounded first derivatives, i.e., ip is locally Lipschitz. He also 
analyzed how ip behaves near the so-called two-particle coalescence points, i.e., 
those points in S(X) where exactly one term in the sums representing V (see 
is unbounded. 

Generalizations with new insights for those points in where more than one 
term in (|1.1|) is singular were obtained in (10, and more recently in [5| . We mention 
that the present authors in [S], and 0] studied the smoothness of the electron 
density, a question related to the present investigation; we shall not discuss this 
further here. 

Suppose we have a solution %p to Hip — Eip, _E e K, with H as in 1)1. 2|) or H1.3|l . 
We want to find a representation for ip) 



(1.3) 



H = Hn{Z) 



A + V 




(1.4) 



HiP ^ EiP, E eR, ipe L2(R3^). 



ip{xi, . . . ,x„) = J^{xi, . .. ,Xn 



) (pixi, . ..,Xn 



such that (p is as smooth as possible and is a universal (i.e., not depending on ip 
or E) positive factor reflecting the behaviour of the potential V near S(X). This 
means that for any two solutions ipi ,ip2 of a fixed Schrodinger operator H1.2|l (or 
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the function IF will be the same, i.e., 

Since it is already known from one-electron atoms that -0 is just locally Lipschitz, 
T cannot be smoother than that. We shall see that by choosing in a special way 
one can say a lot more. Let us first recall some of the ideas developed in P. 

Suppose V' is a solution to (—A + V^)?/; = in C . Set = e^0, then 
satisfies 

(1.5) A0 + 2VF • V0 + (AF + |VFp + (E - Y))<\) = 0. 

Now assume H = — A + is given by 1)1. 2|l . The specific nature of the Coulomb 
potential makes it possible to find an explicit F such that AF = V , namely 

^ N L ^ 
F(x)=F2(x) + - ^ \x,-x,\. 

i=i i=\ i<i<j<N 
We have given F an index 2 to indicate that F2 is a sum of functions each only 
depending on the coordinates of two particles. If we insert F2 into 1)1.5(1 we obtain 

A02 + VF2 • V02 + (|VF2|^ + F)02 - 0, 

where we have also given (j) an index 2 to show that it is associated with F2. The 
regularity properties of (f>2 are now determined by the regularity of V-F2, respec- 
tively, |VF2p. Since VF2 is locally bounded, standard elliptic regularity theory 
(see Sectional) gives us that 

(1.6) 02GCi'"(f7) for aG(0,l). 

(For the definition of the Holder-spaces C'°'", see Definition 12. l|l . Since VF2 is just 
bounded and not continuous, one cannot in general expect anything better than 
(|1.6() . Note that since -0 = e^^02 we have 

(1.7) VV' - (VF2)0 G Cin) for a G (0, 1). 

This is a general formulation of Kato's cusp condition [13| which plays an important 
role in the numerical treatment of (|1.4|) . (Here, and in the sequel, by a 'cusp 
condition' we understand a condition a solution ■0 has to satisfy at a point in the 
singular set E(X)). 

We are now ready to state our main result about the regularity of 0. 

Theorem 1.1. Suppose ip is a solution to Htp = Eijj in C_ M.^^ where H is given 
by ini). Define yu ^ Xi - Xg, i e {1, . . . ,N},e e {1,. ..,L}. Let 

(1.8) F = e-^=+^3 
with 

L N 

(1-9) F2(x) = --^^Z,|y,,| + - 



4 ' 

= 1 4=1 l<i<j<N 



L 

(1.10) F3(x) = Co^ J2 Z,{y,,e■yJ.i)ln{\y^,i\' + \y,J\')), 

1=1 l<i<j<N 



2-7r 
12ir • 

Then 



where Cq 
Then 

(1.11) ^ = ^03 
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with 

(1.12) 03 e c^^\n). 

Furthermore this representation is optimal in the following sense: There is no 
other function T depending only on on X, Z and on N , but not on ijj or E itself 
such that ijj = J-(j) with (j) having more regularity than C^'^{Q,). 

Remark 1.2. 

(i) Of course one can consider more general Hamiltonians, for instance molecular 
Hamiltonians where the nuclei are allowed to move. Kato \V6\ considered this 
case. Our results, suitably m,odified, extend to this situation. We concentrate 
on the model with fixed nuclei since this is the 'standard model ' in molecular 
physics. 

ill) For the proof of Theorem \l.l\ a special regularity result (see Theorem \2.h]) for 
solutions of the Poisson equation Au = g will be vital. Roughly speaking, if 
g € has a certain multiplicative structure, we can show that u g C^'^, and 
not only u € C^'",a £ (0,1) as in general (see Provosition \2.^) . This result 
is of independent interest. 

(iii) Note that each term in the sum F2 is either a term involving the coordinates 
of one electron and one nucleus, or the coordinates of two electrons, whereas 
the terms in F3 involve the coordinates of two electrons and one nucleus. In 
the representation (I1.1U|) of F^ no terms involving the coordinates of three 
electrons occur; see Section\^for details. 

The fact that no terms involving the coordinates of four and more particles 
show up in F^ stems from the fact that in the summands contributing to | Vi^2 P 
only terms involving at most three particle coordinates occur (again, see Sec- 
tion\^for details). 

(iv) An immediate consequence of Theorem is the following sharpening of 

W - V^VFa + VF3) G c"'\n). 

(v) Attempts to approximate many-particle wave functions by a product as in 
()1.11|) are common in computational chemistry and physics. There, such an 
T is usually called a ' J astrow factor'. 

It is also interesting to consider the regularity of ijj near the zero-set Af{ip) = 
{x e I V' = 0} of V'- A simple argument shows that Theorem 11.11 actually 
implies that VV' : ■N'ii') ^ is locally Lipschitz, whereas Vij) is just bounded in 
I](X) \A/'('!/')- By 'locally Lipschitz' we here mean the following: For all closed balls 
K C M=^^, there is a constant C = C{K) such that | V?/'(x) - VV'(y)| < C(i4:)|x-y| 
for all X, y e 7V(V') n K. Indeed, writing Vip = tpy{F2 + F3) + cxp(F2 + F3)\/(j)3, 
we get, for x G Afi-)p), that V'0(x) = exp(i^2(x) + F3(x))V(/)3(x) since V(-F2 + F3) 
is bounded. The assertion follows, since both exp(F2 + F3) and V03 are Lipschitz 
in 

In jSI it was shown for a wide class of potentials that at their zero-sets real valued 
distributional solutions (which for these potentials are then actually continuous 
functions) to (—A + V)u — are, roughly speaking, by one degree smoother than 
away from their zero sets. So the observation above extends these results to the 
Coulombic case. The potentials considered in [S| were of Kato type, K^'^ , where n 
is the dimension (in our case, n = 3A^) and S G (0,2); see ^H] for definitions and 
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many far-reaching results concerning these potentials. In US| (see also JHl) it was 
shown that solutions are locally for 8 < 1 and C^'^''^ for S e (1,2). However, 
since the Coulomb potential V in (|1.1|) is in 

j^3N,s fQj. alU < 1, but not in X^^-i 
these results are not sharp and actually weaker than Kato's result. 

It is not surprising that logarithms occur in H1.10() . Such terms have been consid- 
ered in classical work by Fock 3 for the atomic case; see Morgan JTP for an analysis 
of these 'Fock-expansions' for two-electron atoms. That paper also contains many 
references to earlier work on such expansions. 

Proof of the optimality of the representation Hl.ll|l : It suffices to find a simple 
example. Consider the one electron atom whose Hamiltonian is given on M.^ by 

Z 



H = -A--—r , X = {xi,X2,X3) e 



j>3 



With "01(2;) = e ? 1^1 and 0^2 — xie ? 1^' we have 



"2- 



Write = and ip2 = T(fP-'^ . Now Vi > and if we had an T which would 

allow more regularity of the (^^*^'s, then 

would be better behaved than just C^'^. But near the origin the right hand side 
of H1.13|l behaves like a:i(l -I- ^\x\) and this is just C^'^, i.e., the second derivatives 
are bounded but not continuous. □ 

The results in Theorem 1 1 . 1 1 are not well suited for obtaining a priori estimates. In 
particular neither nor F-j, stay bounded as |x| tends to infinity so that if, say, 
■0 G L^(M^^) then 03 is not necessarily in L^(R'^^). These shortcomings will be 
dealt with below in a similar way as in 1^. 

Definition 1.3. Let x e C^{^), < x < L with 

1 for\x\<l 
for \x\ > 2. 

We define 

(1-15) i^cut — -Pj^cut + -Fb^cut, 

where 

L N 

(1.16) F2,™t(x) = E E \y^^^\ 



(1.14) x{x) 



i=l 1=1 

1 . , 

X ( I -^i Xj\) \Xi Xj\, 



- E 

4 ^ 

l<i<j<N 



(1.17) F3,eut(x) 



L 



Co XI E z^x{\v^A)x{\v3A)iy^,^■yj■^)^^{\y^A^ + \v3A^))^ 

l=\ l<i<j<N 
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and where Co is the constant from (|1.1Q(I . We also introduce (/)3,cut by 
(1.18) V = e-^-'03,cut. 

Theorem 1.4. Suppose ijj is a solution to Hi/j = Ei/j in M.^^ . Then for all < 



R < R' there exists a constant C{R,R'), not depending on ip nor xq G R , such 
that for any second order derivative, 

d^^—^ , z,j = l,2,...,iV, A:,^=l,2,3, 

the following estimate holds: 

(1.19) ||9V-V'5'i^cut||L^(i33«(xo,fl,)) < C(i?,i?')llV^IU-(i33«(xo,fl/)). 

Remark 1.5. Theorem \1.4\ strengthens results obtained in & . More precisely, to 
prove Theorem \1.4\ we will show that 

(1.20) ||</'3,cut||ci.i(-B3iv(xo,fl)) < C'(-R:-R')ll</'3,cut|U~(S3N(xo,-R'))- 

The estimate (I1.19|) is then a trivial consequence of (|1.20|l . ( On the other hand, 
((rT^ and (fL^ imply (fT^ ). 

The estimate H1.20|l is a strengthening of Proposition \1.6l below to a — 1. We 
state and prove the proposition here, since we need it in the proof of l|1.2()|l . R 
essentially follows from ideas in 9 . 

Proposition 1.6. Suppose is a solution to Hi/; = Ei/j in M.^^ . Then for all 
< R < R' and all a € (0, 1) there exists a constant C{a, R, R'), not depending on 
"0 nor Xq G K^^, such that, with 03, cut defined as above, 

(1.21) |i03,cut|!ci.°(B3N(xo,-R)) < C'll'/'3,cut||L°°(S3N(xo,fl'))- 

Proof of Provosition \T~B : Note first that with ip = e^^ ''"' 02, cut, H1.5() and AF2 = V 
gives 

(1.22) A02,cut + 2Vi^2,cut • V02,cut 

+ (A(F2.cut - F2) + |VF2,cut|' + S)02,cut = 0. 

It follows from the form 0/ F2,cut and F2 (see (|l.lt)|l . (I1.14|) . and (|1.9|) ) that the 
coefficients in (|1.22|l above are uniformly bounded in M'^^. Therefore, (|1.21|l . with 
02, cut instead 0/ 03, cut, follows from Provosition \2.'/A To get 11.21|l with 03, cut, note 
that 03, cut — e"^-' ""' 02. cut, o-nd that -Fs.cut G C^'"(]R'^^) and has compact support 
(see (frT7|l and ifLHIl ). □ 

We point out some consequences of Theorem 11.41 which can be viewed as cusp 
conditions for second order derivatives of ^. Indeed, we can relate the singularities 
of the second order derivatives of i^cut with those of the second order derivatives 
of "0 in a precise way, thereby obtaining certain identities. Here we only explicitly 
state some representative cases. 

Corollary 1.7. Let tp be a solution to H'tp = Eip in R^'^ with H given by ^1.2\ . 
(i) Let 1 < i < j < N , and fix any point zq — {zi, . . . , zn) G R'^^ with Zi = Zj = 
z. 

Then 

(1.23) hm (|x,-a;,|V..V,^) + i^(zo) -0. 
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(ii) Let l<i<N,l<£<L, and fix any point zq = {zi, 
Zi = Xi, Zj ^ Xt,i ^ i. 
Then 



,zn) e with 



(1.24) 



lim 

ii-»o 



{\x,~-Xi\ Aiip)+Ze ?/'(zo) 



L~(B3„(zo,i?)) 



= 0. 



Proof: In order to show (|1.23|l we first show that 



(1.25) 



Hm \xi - Xj\ Vi • Vj Fcut(x) 



It suffices to consider the limits for \xi — Xj \ • VjF2 and \xi — Xj \ Vi • VjF^. An 
easy calculation shows that 



lim \xi — Xj\Wi ■ Vj -F2(x) 

x^zo 



li z ^ Xe for all i then V; • Vj F3 is smooth near Zq. We therefore only need to 
consider the case z = Xi. We have 

CaZeV, ■ Vj{{{x, - X^) ■ (x, - Xe)) In {\x, ~ Xef + \x, ~ Xif)) 
= -iCoZi In {\x, - Xif + \xj - X,|2) + 7^, 

where 77 is bounded in a neighbourhood of zq. Noting that 

\x, -x,\<^{\x,- Xtf + \x, - X,|2)'/', 

we see that 

lim V, •VjF3(x) =0. 

x^zo 

Using the triangle inequality we obtain 



\xi - x^\ (Vi • VjV) + ^V'(zo) 



< 



L°°(B3„(zo,-R)) 
- Xj \ ((V, • VjV') - (V, • VjFeut)V') 



Icci - Xj \ (Vj • Vj Fcut)ip + ^-ipi^o) 



This, jnni), and fL^ imply (fT^ . 

The proof of (|1.24() is similar. Just note that 



L~(B3jv(zo,fl)) 

L~(B3jv(zo,-R))' 



□ 
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1.2. Organisation of the paper. For simplicity we shall only give the proofs of 
Theorems 11.11 and 11.41 for the atomic case (i.e., £ — 1, Xi = and Zi = Z, see 
()1.3|) ). Indeed, no additional complications arise for molecules. Also, we only give 
the proof of Theorem II. II in the case il — M.^^ . 

In subsection 1 1 . 31 we define some notation to be used in the entire paper. Sec- 
tion El contains standard elliptic regularity results in subsection 12. II Subsection 12. 21 
contains in particular the elliptic regularity result Theorem 12. 61 which is proved in 
subsection l2.3l Theorem l2.6l is the essential new mathematical input necessary for 
the proofs of Theorems 11.11 and 1 1 . 41 These proofs are given in Section 13 — the proof 
of Theorem 11.11 in subsection 13.11 and that of Theorem 11.41 in subsection 13.21 The 
Appendices EI ^ and |0 contain the construction of solutions to certain Poisson 
equations. These solutions is another important ingredient for the proofs of the 
main theorems. 



1.3. Notation. Throughout the paper, constants occuring in inequalities will be 
denoted by the symbol C, although their actual value might change from line to 
line. 

For X G M" {n > 2) we write x = rui, with r — \x\, lo — x/\x\ € S""""", the unit 
sphere in M". Denote by Bn{x, r) the open ball of radius r > around x. 

We denote by Yi^m(ti;) the normalised (in L^(S"~^)) real valued spherical har- 
monics of degree I, I G Nq, with to = 1, ... , h{l) — 1, where 

^^^^ ^ ■ 

Then {Yi^rn}i,m constitutes an orthonormal basis in i^(§"^^). 

The Yj^m's are the eigenfunctions for C^, the Laplace-Beltrami operator on §"^^: 

where — ^ is the angular part of the Laplacian in M" , so 

n-1 d 



- A = - 



dr^ r dr 

We define 'P/^ to be the orthogonal projection in L^(§"^^) on Yi.m,: 
and 



(1.27) rl-^ = Yl ^S- 

We denote f)|"^ = Ran(p/"^). 

By abuse of notation, for a function / : K" ^ C we write f{ruj) = f{x), and, 
whenever /(rp •) : S"^^ ^ C is in L^(§"~-^) for some Vq G (0, oo), we write 

(^/m/)(^ot^) = Yi^miuj) / y;,m(w)/(row) duj = fi,m{ra)Yi_jn{uj). 
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2. Elliptic regularity 

In this section we collect results on regularity of solutions to second order elliptic 
equations needed for the proof of Theorems 11.11 and 11.41 The results fall in two 
parts, known ones (in subsection I2.1|l and new ones, developed for our purpose, 
and of interest in themselves. The latter are in subsection 12.21 The result of main 
interest is Theorem 12.61 which is proved in subsection 12. 31 

2.1. Known results. 

We start by recalling the definition of Holder continuity: 

Definition 2.1. Let ft be a domain in M", fc e N, and a € (0,1]. We say that 

a function u belongs to C'°'"(0) whenever u e C^[Q), and for all (3 S N" with 
\(3\ = k, and all open balls Bn{xo,r) with Bn{xo,r) C fl, we have 

iDf^uix) - D'^u(y)\ ^, 
sup J _ |, < C{x^.r). 

For any domain D,' , with fl' CC il, we define the following norms: 

l/3|<fc 

\Dl^u(x) - D/^uiyM 

[u]k,a,n' = _ Ic. ■ 

l^l^j, a:,yen', \x y\ 

For k = we use the notation C°'{fl) = C'^'"(ri) and [u]a,n' = [u]a,a,n'- 
Furthermore, for a function u g C"(M" \ {0}) we define 

(2.1) Ikllcocs"-!) = sup |m| + M„_s„-i, 

\u{x) - u{y)\ 
= sup , .o. ■ 

x,y&"-^,x=ty 1-^ y\ 

We will need the following result on elliptic regularity in order to conclude that 
the solutions of elliptic second order equations with bounded coefficients are C^'". 
The proposition is a reformulation of Corollary 8.36 in Gilbarg and Trudinger [Jj, 
adapted for our purposes: 

Proposition 2.2. Let f^o be a bounded domain in M" and suppose u € W^''^{VIq) 
is a weak solution of Au + ^j^j'^ + T^"" — 9 in flo, where bj, W,gG i°°(ilo). 

Then u £ C^'°'{VIq) for all a G (0, 1) and for any domains £7', il, fJ' C fi, C fio 
we have 

\\u\\c^,<^(n') < C(sup|u| +sup|g|) 
n SI 

for C = C{a, n, M, dist(17', dfl)), with 

max{l, max l|W^llL~(n), I15||l-(0)} < M. 

J — i , . . . , Tl 

We further need results concerning the regularity of solutions of the Poisson 
equation. These regularity properties are based on the regularity properties of the 
Newton potential of the inhomogeneity. For our further considerations we recall 
here the properties of this function. 
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Let g g L°°(ri) for H. a bounded domain in M", n > 2. The Newton potential of 
g is the function w defined on M" by 

(2.2) wix) = [ T{x-y)g{y)dy 

with 




From 15, Theorem 10.2 and 10.3] we have 

Proposition 2.3. Let Vl C W\n >2, be a bounded domain, then: 

(i) Ifg€ L°°{fi), then w G C'^'"{n) for all a G (0, 1), and Aw ^ g in n holds in 
the distributional sense. 

(ii) Ifge C'='"(f]) for some keN and some a G (0, 1), then w G C''+^'"{n). 

Since every solution to the Poisson equation can be written as a sum of the 
Newton potential of the inhomogeneity and a harmonic function, the above implies 
in particular the following well-known result: 

Proposition 2.4. Let g G C'=^"(f^o) for some fc G N and some a G (0,1), and 
assume u is a weak solution to Au — g in D,q. 

Then u G C'^+^'"(Slo). Furthermore, for any domains ft', ft, fl' (Z il, CZ ilo, 
we have 

(2.3) ||u||cfc+2.o(j^,) < C(sup|u| + ||g||c'=.°(n)), 

n 

with C = C{n, k, a, dist(f7', dQ)). 

The next lemma, which is taken from Gilbarg and Trudinger T, Lemma 4.2], is 
essential for the proof of the main regularity result in subsection 12.21 

Lemma 2.5. Let be a bounded domain in M", n> 2 and let g G C°'{fl) r)L°°{il) 
for some a G (0, 1] . 

Then the Newton potential w of g (given in (|2.2I) ) satifies, for x d Q and i,j = 
1,2, ...,n, 

D,jw{x) = / D,jT{x - y){giy) - gix)) dy 
J no 

(2.4) ~gi^) f D,T{x~y)v,{y)da{y). 

Jdno 

Here, fio *s a*^?/ bounded domain containing for which the divergence theorem 
holds, and g is extended to vanish outside fl. In the last integral, da denotes the 
surface measure ofdfla, and Vj the j-th coordinate of its (outwards directed) normal 
vector. 

2.2. New results. 

We here collect a number of more explicit regularity results needed in the proof 
of Theorems 11.11 and 11.41 

The following result shows that one can push the C^'",q; G (0, 1), in Proposi- 
tion l2.2l to C^'^ in certain cases. 



SHARP REGULARITY FOR WAVE FUNCTIONS 



11 



Theorem 2.6. Let g E L°°(M'^), k > 2, be a homogeneous function of degree 
which has the properties g G C"(M'' \ {0}) and g\sk-i is orthogonal to l)^'"' (the 
subspace of L^{E>'^~^) spanned by the spherical harmonics of degree 2). Let f S 
C°'(U.'^) for some d> and let u G C^'°'{MJ'~^'^) be a weak solution of the equation 

(2.5) Au{x',x")^g{x')f{x") 

where x' eR'', x" € M'^, A = A^- + A^// . 

Then u G Wf^^ {MP'), n = k + d, and the following a priori estimate holds: 
For all balls Bn{z, R) and B„iz, Ri) in M" where < R < Ri, z E W\ 

!c°(7rdB„(z,i?i)) 



sup l-Dy uj < C ( sup \u\ + ( sup \g\ ) 



(2.6) +( sup I/I) ||g||c»(s-i)) 

7TaB^(z,Ri) ' 

With C = C(n, a, R, Ri). Here naix' , x") = x" for x' E R'', x" E for d > 0; for 
d = 0,7rrf(a;') = 0. 

Remark 2.7. 

(i) The case d — means that f is a constant and the terms in (|2.6I) with f then 
equal this constant. 

(ii) The reason for the condition k > 2 will become clear in the proof of the 
theorem, when Lemma \2. 51 its applied. 

(iii) Note that if k = Q,d > 2, one has stronger conclusions: Equation (|2.5|l be- 
comes Auiy) = f{y) with f E C"(M''), so by Proposition^^ u E C'^^°'{W^). 
The a priori estimate analogous to H2.6|l is then a consequence of Holder- 
estimates for u (see e. g., [7| Corollary 6.3]/ 

(iv) Using the standard fact (H Theorem 4 in 5.8]; that Wi^'~{R") = Ci^o'c (K") 
(with equivalent norms) we may replace the ierm sup^^^^ b) \Diju\ by [u]! i s„(2,i?) 
071 the left hand side in (|2.6I) . 

(v) For the special solution to 12.51) given by the Newton potential of gf, the 
estimate H2.6|l holds without the term sup^^^^^^^j |u| on the right hand side 
(see (j^TOjl ). 

Since the proof of Theorem 12.61 is a bit lenghty we present it separately in sub- 
section |^1 

The following proposition, on solutions to Poisson's equation, when the inhomo- 
geneity / in Au = / is a homogeneous function, is needed often in the paper. 

Proposition 2.8. Assume that the function g satisfies 
g{ru) = r^G{uj) with G E (§""!) and vi%G = 0. 
Then there exists a solution u to 

(2.7) Au = .g on B„(0,i?)cM", 
satisfying u{ruo) = r^+'^U{uj) with U E Ci'"(§"-i) for all a E (0, 1). 

Proof: Let 

gi,m{r) ^ / g{ruj)Yi.^{uj)duj = r'^ I G{uj)Yi.^{uj)duj = r''gi^rn- 
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Then (see (fT^ for h{l)) 

oo h{l)-l 
l=0,l^k+2 m=0 

since gk+2,m = for all m. 
Now define 

oo Hl)-1 
l=0,l^k+2 m=0 ' ^ 

with k) = {k + 2)((fc + 2) + n - 2) - + n - 2). Note that bi{n, k)) ^ for 
Z ^ /c + 2. Since e i2(§«-i) (gi^(,g ^ i°°(S"-i)) the sum ^ 

therefore converges in L^(S"^^). 

Make the 'Ansatz' u{ruj) = r^+'^U{uj), and denote for e N 

N h(l)-l 
gjy{ruj) ^ ^ ^ 9Lmr^Yi^rn{^), 
i=0,i^fe+2 m=0 
AT 

l=0d^k+2 m=0 ' ^ 

Now let e C^{B„{0, i?)), then, using that C^Yi^rn ^ l{l + n ~ 2)11,™, 

(2.9) / (j){Au — g) dx = / (A0)(u — u^v) c^a; + / (j)[g n — g) dx . 
Jb„{q,b.) JB„(o,i?) Jb„{o,b.) 

Since u — ua? — > 0, g — ^ (in - sense) for N ^ 0, the RHS of 12. 91 tends to 
zero for 0. Hence u = r^'^'^U{u}) solves IT7I in the distributional sense. With 

w the Newton potential corresponding to g fsee l2.2f) . we have w G C^'"(i?„(0, i?)) 
due to Proposition l2.3l and u — w is harmonic, so u € C^'"(i?„(0, R)). This implies 
that [/ e Ci'"(§"-i). □ 

We prove the following useful lemma: 

Lemma 2.9. Let G : U ^ M" for U C M"+™ a neighbourhood of a point (0, j/o) e 
M" X «"\ Assume G(0, y) = /or all y such that (0, y) G [/. Lei 

^•G(x,2/) x^O, 
x = 0. 



Then, for a G (0, 1], 

(2.10) G e C^-"{U] M") / e C°'"(t/). 

Furthermore, ||/||c"(c/) < 2||G||c°((7)- 

Proo/: Let a € (0, 1]. We need to estimate f/^^'^l^~/^'"'''^f2 . 

■' V ' J |(a:i,i/l)-(2:2, 1/2)1° 

Suppose first that X2 = 0. Then f(x2, 2/2) = and we get 



|/(xi,2;i)-/(0,y2)| ^ 



F^-G(xi,yi) 



|(a;i,yi)-(0,y2)h " kil" 
< l|G||c"(c/), 
since G € G"(t/;R") and G(0,2/i) = 0. 



< 



Xl 



\G{x^,yi)\ 
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Next, suppose < 1 < \xi\. By the triangle inequality: 

I "^1 

\f{xi,yi) - f{x2,y2)\ < 1^ • {G{xi,yi) - G{x2,y2)) 



+ 



Xl 



X2 



) ■ G(x2,J/2) 



\Xi\ \X2\ 

Using that G is C" and that G{0, 1/2) = 0, we get 

|/(a;i,yi) - f{x2,y2)\ 

< WGWciu) ilixuyi) - (a;2,y2)r + {-p-r - 7^ 

' ' y ^\xi\ \x2\ 

To control the last term — divided by \{xi,yi) — (a;2,y2)|" — we first derive a lower 
bound on \{xi,yi) ~ {x2,y2)r- 



- (a;2,y2)| > \xi ~ X2\" 



X2 



= {\X,\ - \X2\)' + \X,\ |X2|(^ - ^) ^ l^ll 

Therefore, using the assumption < 1x21 < 

\ixi,yi) ~ (0:2,^2)1" > 1x21" 
This finishes the proof of the lemma. 

The following obvious lemma is used repeatedly throughout the paper: 



Xl 



X2 
1^2 I 



□ 



Lemma 2.10. Assume f{rLu) ^ r^G{Lu) with G e G'^'^iW \{0})n L°°{W). Then 

f e Ci'i(R"). 

Proof: The first derivatives of / trivially exist and are continuous. Therefore it 
suffices to show that all derivatives of / of second order belong to Lj^^(M"); the 
result then follows from Remark 12.71 lpv)l . 



dxjdx 



or ^ ^G dG 

- = 2<5,,,GH + 2(x,— +x.— 



, d'^G 

dxjdxk 



since G € C^'\W \ {0}). This proves the lemma. 



□ 



Note that better regularity cannot be expected without assuming continuity of G 
at X = 0. On the other hand, if G only depends on a; S and G is continuous 

at X = 0, then G is a constant. 

2.3. Proof of Theorem 1231 

We first investigate, for xq € Bn{z,Ri), the behaviour of the Newton potential 
w as given in H2.2|l . namely 



(2.11) 



w{xo) 



T{xo^y)g{y')f{y")dy 



with y = (y',y") G M''+'* 
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Since u and w are C^'" - solutions of (|2.5|l in Bn{z,Ri) (see Proposition 12. 2|l . 
h = u — w is harmonic. Any harmonic function h in a bounded domain satisfies 
the foUowing a priori estimate (see Theorem 2.10]): 

(2.12) sup|Aj/i| < ^sup|/i| , i,je {l,...,n}, 

K On 

with K compact, C f7 C M", and (5 = dist(i4:, dil). So, by l|TTT|l and ^A^, for 

a;o G Bn{z, R) (recall that h — u ~ w) 

C{n) 

{Ri ~ RY ^b„{zMi) 



|AjU(xo)| < — 2( sup 



+ C(n,i?i)(sup|.g|)( sup |/|) 

(2.13) +\D,jw(xo)\. 

Therefore to prove that u G W^^^{M") and that u satisfies H2.6|l it obviously 
remains to show that w satisfies the a priori estimate (|2.6|l . This will be done via 
Lemma 12.51 and will finish the proof of Theorem 12.61 

We proceed as follows: Define N = {{x',x") G R" \ x' — 0} and note that 
|iV| = (|A^| denotes n-dimensional Lebesgue measure of N) and that for every 
ball Bn C M", Bn \N is still a domain. For this the assumption fc > 2 is vital (see 
also Remark i2. 71 (|u|l). Note also that (still with xq E Bn{z,Ri)) w can be written 
as 

(2.14) w{xo)^ f T{xo-y)g{y')f{y")dy. 

Jb„(z,Ri}\N 

Taking into account the Holder continuity assumptions on g and / it is easily seen 
that for every domain O C M", gf G C°(f} \ N). Hence (|TT^ and Proposition ITl 
implies that w G C^^" (S„(z, i?i) \ TV) . 

Now we are ready to apply Lemma 12.51 Pick fl ~ Bn{z,Ri) \ N and JIq = 
Bn{z, R2) with Ri < i?2, then we obtain from (|2.4|l . for xq G Bn{z, Ri) \ N, that 

D^jw{xo) = / D,jT{xQ - y){{gf){y) - igf){xQ)) dy 

JB„{z,R2) 

- {9f){xa) I DiT{xa ~ y)^j{y) d(T{y) 

JdB^(z,R2) 

(2.15) =I{xo) + Jixo). 

Here as before gf is extended by zero outside Bn{z,Ri) \ N. Noting again that 
|A^| = 0, we can use this integral representation to derive the a priori estimates on 
DijW. We want to show that for < i? < i?i 

(2.16) sup \Dijw\ 

B„iz,R) 



< C2 



( sup \g\) ||/||c°(7r,B„(2,i?,i)) + ( sup I/I) ||5llc°(S'=-i) 

S*--! TTdB^(z,Rl) 



where C2 = C2(n, a, i?i — R,R/Ri). Inequality (|2.16(l together with inequality 
(|2.13(l will yield the desired a priori estimate (|2.6(l and implies in particular that u G 
M^j^'^°°(R'^+'^). So to finish the proof of Theorem 12. 61 it remains to prove inequality 
(|2.16|l . For this we have to estimate the integrals I{xo) and J{xo) in l|2.15|l . We state 
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the estimates in the following lemma fLemnia l2.11|l . which we then apply to prove 
inequality (|2.16|l . The proof of Lemma l2.11l is given afterwards. For convenience we 
shall henceforth use the following notation: B = Bn{z, R), Bj ~ Bn{z, Rj), j = 1, 2. 

Lemma 2.11. With /(xq) and J{xo) as in H2.15|) we have the estimates 
\I{xo)\ < C{n) sup \g\) ( sup |/|) 

+ C{n,a) {R, ^ Rr{ sup \g\)[f]^ 
+ C{n,a) ( sup I/I) ||g||c°(S'=-i). 

R2 



(2.17) 
(2.18) 



|J(xo)|<C(n) sup \gf\ 

B„{z,R) 



R2 — R 



n-l 



for xo e B„{z,R)\N. 



Combining the inequalities 1)2.17(1 and ((2.18(1 with ((2.15(1 leads to the a priori 
estimate 



(2.19) sup \D,jw{xo)\ 

xoeBn{z,R) 

R 



+ C(»,a)f(fl, -if)°(sup|s|) [/]<,,,B +(»P I/I) llfllc-cs.-.) 



TTdBi 



Finally we pick R2 — 2Ri and obtain, with C — C{n, a, i?, i?i), 



(2.20) sup \D,jw\ < C 

B„{z,R) 



. sup I5I) ||/||c°(^<iSi) 



, sup I/I) ||5llc°(s'=-i) 



This finishes the proof of 1(2.16(1 and according to our previous considerations the 
proof of Theorem 12. 61 It remains to prove the estimates in Lemma [2. Ill 

Proof of Lemma In. 111 . We start by proving the estimate 1(2.18(1 on J{xo). For 
y G dB2 and xq e Bn(z, R)\ N we have |a;o — y\ > R2 — R- This, and 

C{n) 



\D,Tixo~y)\ < 



Fo - y 



n+l ■ 



yields 



\JM\ < \igf)ixo)\ 



|Ar(a:o - y) ujiy) da{y)\ 



dBo 



<C{n) sup \gf\ 

B„{z,R) 



i?2 



i?2 — R 



verifying ((2.1811 . 

It remains to prove the estimate ((2.17(1 on I{xo). This is more involved. With 

R' = Ri- R a.ndn ^ Bn{z, R2) \ N, write 

(2.21) I{xo)=h{xo,R')+l2ixo,R') 
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with 

/i(xo,i?') = / A,r(xo - v){{9f)[y) - (.g/)(a;o)) dy, 

Jn\B„(xo,B') 

/2(xo,i?')= / D,,T{xa-v){{9r){v)-{9f){xa))dv. 

Jb„(xoM') 

Clearly we have (since g/ = on ^2 \ Bi and g is homogeneous) 
|/i(^o,i?')l <2(sup|5/|) / \D,,T{xo-y)\dy 

Bi Jn\B„(xo,R') 

(2.22) <C(n)(sup|g|) (sup|/|)(§)". 

The estimate for l2{xo, R') will be more involved and we need several steps. 
First notice that 

(2-23) (A,r)(x) ^ ^ , 

where P2 is a homogeneous harmonic polynomial of degree 2 (which clearly depends 
on the indices we suppress these for simplicity). Use polar coordinates x — 
ruj, r — \x\, uo — x/\x\, and obtain (using /g„-i P2{oj) duj — 0, and l|2.23|l ) that 

(2.24) h{xo,R!)=C{n) ( [ r'^ P2{u;){gf){xo + ru;) du; dr. 

Jo JS"-i 

Denote x €W hy x = (x', a;") = nu = r(w',a;") where G M'', cj" e R'^ (so that 
|w'|2 _)_ |t^"|2 = x- ^hen d = 0, w" = 0). With this, write 

(2.25) hixo, R')=C{n) ^ r-\l^^\xo,r) + I^^\x^,r)) dr 

Jo 

with (4^^ = when d = 0) 

4\xo,r) = f P2{cj)g{x'o + rcj'){f{x'i + ru;")-f{x'^))du;, 

lP{xo,r) = /«) / P2{L^)g{x'o + rLu')dLu. 



We need to estimate |/2^''| and 1/2^'! such that we gain a suitable r-behaviour for 
small, respectively, large r which will enable us to estimate Ihi^o, R')\- 

Firstly, due to Lemma 12.51 {gf){y) is defined to be zero for y G B2 \ Bi in 
I{xo) and formula H2.15|l holds for xq G Bi \ N. Using this formula just for xq G 
Bn{z, R)\N we have xo + roj e Bi for all r with 0<r<R' = Ri—R and therefore 
(up to the zero set N U dBi) we can make use of the Holder continuity properties 
of g and / for the points xq and xq + ruj in the integrals above. 

Using the Holder continuity of / and the homogeneity of g we obtain 



l4"ta,r)| 



S"-i r \uj I 

< C(n)r"( sup I5I) [f]c,n^Bi 
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(for [f]a,TTaBi, see Definition 12. 1(1 . Hence 
(2.26) 



r-i|/«(xo,r)|dr < ^ (i?')"( sup \g\) [fU^.B,. 



To estimate the second term in (|2.25|l . we write in the fohowing Xq — \xQ\r] with 
rj € S*"'^^ and define s by r = |a;o|s. Then 

(2.27) 



r "'"/2^''(a;o, r) dr 



1/(4) 



R' 



P2('J-') (7(|xq|77 + ruj') duj dr 



P2{lo) gii] + suj') duj ds 



duj ds. 



where we used that g is homogeneous of degree zero and \x'^^\ ^ 0. Because of 
the s~^-term in the s-integral we have to control the cj-integral for s — > and for 

Define, for < si < S2 < oo, 

(2.28) K{si,S2)= f\-' f P2{uj)g{f] + suj') 

The behaviour of K for different regimes of si and S2 is expressed in Lemma l2.12l be- 
low. Applying it, we get that (for all |a;o| e (0, cxd)) |/ir(0, i?'/|xQ|) I <C{n,a) ||g||c°(s'«-i) 
(for ||g||c°(s''-i): see Definition 12. 
Since due to (jTTrjl 

r r~^lP{xo,r)dr = \f{x'^)\ \K{0, R' /\x'o\)\ 
Jo 

we obtain 
(2.29) 

Further via l|223l, and itT^ lead to 



R' 



r i/f^(xo,r)dr < C(n, a) ( sup |/|) ||5llc"(S'=-i)- 



(2.30) 



ixo,R')\<Cin,a) (i?')" ( sup I5I) [/]„ ^^^^ 
+ C{n,a) ( sup I/I) bllc^is"-!)- 



The estimate (|?T7ll now follows from and (|On|) . 

Proving Lemma [2.121 below will finish the proof of Lemma [2. Ill 

Lemma 2.12. With K as in (|2.28|l we have: 
(i) si = 0,S2 < 1/2.- 



□ 



(2.31) 

(ii) 1/2 < si < S2 < 4; 
(2.32) 



|-f'^(0,S2)| < C(n, Q;)[g]„_s)=-i. 



^^(51,52) < C(n) sup |5|. 
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(iii) 4 < Si and S2 oo: 



(2.33) 



|^(si,S2)| < C(n,a)[5]^gfc_i + C{n) sup 



Proof of Lemma l2.1Sl : 

(jij: Since /g„-i P2{i^) duj = and g is homogeneous of degree we have 



if(0,S2) 



Note that \ri + slu'\ > 1 - s > 1/2. Since g e C"(S'="i) we obtain 



|if(0,S2)| <C(n) [5] 





rj + suj' 


Jo 


7y + sw' ^ 



This, and 



77 + suj' 



< 



s + |l — |r/ + sw'l 
I77 + suj'l 



< 



2s 



ds. 



2s 

< < 4s 



\ri + su}'\ 1 — s 



I?? + suj'\ 
impUes (|2.31|l . 

This follows directly from the definition of K (see (|2.28|) '). 
(|mj: This is the most involved case. We write the unit sphere §"^^ as the union 

of 



(2.34) 



1 



E(s) = {ujc S*"-! I s\J\ < VJ} - {cj e S*"-! I |cj'| < ^ } 



and its complement S(s)'^ (when d — {), I](s) = for s > 1) and write if(si, S2) — 
Ai+ A2 + A3 where (when d = 0, ^2 = for si > 1) 

Ai^ I s^^ ( P2{uj) {g(ri + suj') - g{suj'))dujds, 

Jsi 



A2 
A3 



P2{uj) {g{rj + suj') — g{suj')) dujds, 



S(s) 



P2{uj) g{suj') dhj ds. 



The estimate (|2.33ll is a direct consequence of the following lemma. Proving it will 
finish the proof of Lemma [2. 121 □ 



Lemma 2.13. We have 

(2.35) 
(2.36) 

(2.37) 



1^1 1 < C(n,a)[g]^ 

|^2| <C(n) sup Iff I, 
gfc-i 

A3 = 0. 



Proof: Ai : Note first that since s\uj'\ > 2 and jry + suj'\ > 1 in S(s)'^ we obtain, 
using the homogeneity of degree zero of g and the Holder continuity of g on , 
that 



l^il <C(n)[5]^ ,,_, 



rj + suj' suj' 



\rj + suj'\ \suj'\ 



duj ds. 
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Then by using the triangle inequahty and that s\uj'\ > y/s > 2, we get 
r/ + soj' soj' 



\rj + suj' 



^1+1 \suj'\ - \-q + sijj'\\ 
~ |77 + sa;'| 

< ^ < ^ < ^ 
~ \i-l + suj'\ ~ \/s-l ~ \/s 



which leads to 
verifying (|2.35(l . 

A2 : For d — A2 — Q. For d > 0, the estimate (|2.36|l is a consequence of the 
following lemma, which is not hard to prove using polar coordinates in M" (we omit 
the proof): 

Lemma 2.14. Let |S(s)| denote the n— 1-dimensional surface measure o/S(s). 
Then 

(2.38) |S(s)| < C(n)s-i/^ 
From H2.38|l we immediately get H2.36f) : 

1^2! < f / s-MS(s)| ds) Cin) sup \g\ < C{n) sup 
VJ4 / §fc-i gfc-i 

A3 : We have 

(2.39) A3 = 

as a consequence of the lemma below (when d — 0, (|2.39() is trivially true, due 
to the assumptions on g), since, by assumption, .g|sfc-i is orthogonal to 1)2'^' (the 
subspace of L^(S'^^^) spanned by the spherical harmonics of degree 2). 

Lemma 2.15. Let < k <n and suppose 4> G L'^(S*''^^) is orthogonal to l^j'^''. 
Let 4> denote the following 'natural' extension of 4>: 

^\{x,y)\) \q for\x\^Q. 
Then (j) G L'^i^'^^^) and is orthogonal to ^^2^- 

Proof: Since </> can be expanded in the natural basis of L^(S'^~^) it suffices to con- 
sider a 4) which is the restriction to §'''~^ of a harmonic, homogeneous polynomial 
Ps of degree s ^ 2. Then Ps{x,y) — Ps{x) for {x,y) G is a harmonic homoge- 
neous polynomial in R" of degree s 7^ 2. Therefore 4>, being the restriction of Pg to 
§"-\ is orthogonal in (§«-!) to ()^"\ □ 

This finishes the proof of Lemma 12.131 and therefore finally the proof of Theo- 
rem El □ 



3. Proofs of Theorems II . II and 11.41 



We recall that for notational simplicity we shall give the proofs of Theorcms ll.il 
and ll.4l onlv for the atomic case. 
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3.1. Proof of Theorem ITTTl Let if; satisfy {H - V)i' = in R^'^ , with V as in 
Hl.l|l . and let F2 and F3 be given as in l|1.9|l and Hl.l()|l . Define 03 by the equation 
ip = e^^~^^^(j)3. Recall that AF2 = V. We now make use of Lemma lb.ll below which, 
together with Theoreni l2.6l is the main ingredient in the proof of Theorem ll.il Due 
to this lemma, there exists a function K3 : M."^^ — > R such that AK3 — — |Vi^2p, 
and G3 = K3-F3e C^^^{R^^). Define Ca by 

(3.1) ^ = e^^+^^Cs 

that is, Cs — e '~^^4'3- Since G3 € C^'^(M^^), it remains to prove Lemma ISTTl below 
and that Cs e C^^^{R^^), then 03 e Ci^^CM^^^) wiU follow. 

Lemma 3.1. There exists a function G3 : R^^ ^ R, G3 e C^'\R^^) such that 
the function 

^3(x) ^ Kaixi, . . .,xn) = Z ^"^ {xj ■ Xk)\n{x'^, + xl) 

iZTT ^ — ' 

l<j<k<N 

(3.2) +G3(x) 

solves the equation AK3 = — |Vi^2p; with F2 as in (|1.9() . 

Remark 3.2. Note that the function (x ■ y) ln(x^ + y"^) belongs to C^'"(R^) for all 



N N-l 



a G (0,1), hut not to C^^^' 
Proof: Note that 

(-^-^^ \7F---(^ \ I (sr^ xi-xj XN ~ xj \ 

^ ' ' 2l|xi|"--'|x^|J""4l^|xi-x,|"--'^ \xN-x,\)' 

so that 

,^r.,2 fNZ^ N{N-1)\ Z V- / 

i<j<k<N 



I j3{Xj, Xk,Xl) 



l<j<k<l<N 

(3.4) =ri + r2(x) + r3(x) 



X y \ X — y 



x\ \x-y\ 

X — y X ^ z y ^ X y 



with (x, y,z £ 
(3.5) l2{x,y) 

\x-y\ \x-z\ ' \y-x\ \y - z\ ' \z - x\ \z - y\' 
Therefore it is natural to make the 'Ansatz' 

K3 = fi + k + 0, 

and look for fi, k,h' solving 

Afi^-Ti , AK = -r2 , Aj> = -r3. 
First, it is easily seen that with /i(a;) = G R^, the function 

1 / ^ ^2 1 

j=l l<j<k<N 
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satisfies A/t = -Fi,/! e C°^(M-''^). 

Further, it suffices to find functions k and v sucli tliat 

y) = ■ y) log(a:' + y') + ni{x, y) , Kie C7i'i(R6), 

OTT 

(3.6) with (Aj; + Ay)K(a;,?/) = 72(x,2/), 
and 1/ e C1'1(K9) with 

(3.7) {A^ + \ + ^z)'^{x, y, z) = 73(3;, y, z), 
since letting 

= J X! '^(^i'^;*;) , i>(x) = -^ ^ i/(a;j,Xfc,a;;) 

\<j<k<.N l<j<k<l<N 

gives (A = Ef=iA,) 

^'^W^I Y ((Aj+Afc)K)(xj,Xfc) = -r2(x), 

l<J</£<Af 

l<j<k<l<N 

The functions k and are constructed in Appendices ^ and ^ Lemma 13.11 then 
follows from Lemma FA. II and Lemma FB. II □ 

Remark 3.3. Summarizing, one can say that only those points where the coor- 
dinates of (at least) 2 electrons coincide with that of the nucleus (xi — Xj = 0) 
give rise to the logarithmic terms in K^. These terms stem from the function n 
and are due to the type of singularity of the j2-terms in |Vi^2p- There is no such 
contribution from the function v, i.e., from the j^-terms in |VJ2p- This is due to 
the permutational symmetry of v with respect to the electron coordinates as will be 
seen from the proof of Lemma \B.l\ 

To finish the proof of Theorem II. II it remains to prove that ^3 e C^'^(]R'^^). 
Using {H ~ E)ip = and iJ = — A + V^, we get the following equation for C3 (see 
(tT31) and inU; set F ^ F2 + K3 and = Ca) 

(3.8) AC3 + 2V(F2+i^3) • VC3 

+ (^A{F2+K3) + \V{F2+K3)\^ + iE-V)y3 = 0. 
Using AF2 — V and A_ft'3 = — IVF2P, this reduces to the equation 

(3.9) AC3 + 2V(F2+i^3) • VC3 

+ (|Vif3|' + 2V^^2 • Vif3 +E)C3^ 0. 

This eliminated one of the terms in the equation for ^3 that was only in i°°(]R^^), 
and not continuous, namely |VF2p. 

To deal with the two remaining ones (containing Vi^2), re-arrange the equa- 
tion (|^ : 

(3.10) AC3 + VF2 • (2VC3 + 2C3Vif3) 

+ (IV/C3I' + + ■ VC3 = 0. 
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Define * = (^-i, . . . , ^-jv) : K^^ M^^ by 

(3.11) =2VC3 + 2C3Vi^3- 



That is, = (^'jm, ^'j,2, ^'i.a) : ~* witli 

(3.12) , = 2^^ + 2(,^ , J e {1, . . . , iV} , z € {1, 2, 3}. 
Tiien 

(3.13) VF2 • (2VC3 + 2C3Vi^3) = ViJ^2 • 

Since ^^3,(3 e C1'"(R3^') for all a e (0,1), we have G C"(M3JV) f^^. ^^jj ■ ^ 
{1, . . . , TV}, i e {1, 2, 3} and a e (0, 1). 

Next, let i : R3(Ar-i) ^ M be defined by 

(3.14) *j,i(a;i, • . . ,a;j_i,Xj+i, . . . ,a;Ar) = *j,i(xi, . . . , x^.i, 0, x^+i, . . .,xn), 
that is, by setting Xj equal to zero in 



Furthermore, define, for j < k, j, k G {1, . . . , N}, the functions (f>( - ■ • 



by 

(3.15) <I>(j_fc)(xi, . . . ,a;Ar) = 

^j{xi, . . . ,Xj-i, \{xj + Xk),Xj+i, . . . ,Xk-i, \{xj + Xk),Xk+i,- ■ ■ ,Xn) 

^k{xi, . . . , Xj-i, \ {xj + x/c),Xj+i, . . . , Xk-i, \ {xj + Xfc), a;fc+i, . .. ,xn)- 
The proof of Theorem 1 1 . 1 1 will follow from the following two lemmas: 

Lemma 3.4. Let j,i^^(jk),i be defined according to (|3.12|l . (|3.14|l and (|3.15|l . 

j4sSMme the functions Uj,i,V(^j k)^i solve the equations 

(3.16) Az.,,, = VP,-,, 

(3-17) Ai;(,-fe),, = - 4^-3^ ^U.k),^■ 

Then G C1'1(M3A^). 

Lemma 3.5. Let 'i>j,^(^j k) be defined according to H3.12|l . H3.14|l and H3.15|l . Then 
the functions 

1 Xj Xk 

4 - Xk\ 
Z xj_ 

all belong to C°(M^^) for all a G (0, 1). 

Let us first finish the proof of Theorem ll.il using the two lemmas. 
Let the function U : M^^' ^ M be defined by 

3 w 3 

(3-20) c^ = EE"^-^*+E E 

4=1 j = l 1=1 l<j<k<N 



(3-19) f S ■ (*) - *i) 
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with the functions Uj.i,V(j^k),i solving the equations (|3.16(l and l|3.17(l . Then 

and, due to Lemma 1X1 U e C^' ^R^^ ). 

Let = C3 - f/, then due to (|XTU|l . (E^U, and the form of VF2 (see (jX^ l 

= E - 2"ll ■ " "^^^ " ('^^^'' + ^)C3 - 2Vi^3 • VC3 

Using the fact that X3,C3 G C1'"(M3^), and Lemma 1X31 we conclude that the 
RHS in belongs to C"(M3^) for aU a e (0, 1) . Due to Proposition[0 W <E 

C2'"(]R3W) for aUa e (0,1), and so Cs + e C^^\R^^) (since U e C1'1(M3A^) 
as mentioned above). 

This finishes the proof that C3 e (^^'^(M^^), and therefore ^3 = e'^^C3 e 
Ci'i(R3^), since G3 e C1'1(M3W). 

To finish the proof of Theorem ll.il it therefore remains to prove Lemma m and 
Lemma 13.51 

Proof of Lemma \3.4\ - Firstly, for Uj,i, this is a straightforward application of The- 
oremESl with fc = 3, d = 3{N - 1) and 



X 



g = T^^ , x' = x,e M^ 

/ = |*j. , x" = (a:i,...,x,_i,x,+i,...,a;Ar)eM3(^-i). 

It has already been noted that 'i'j^i E C°'{M.^^) for all a G (0,1) and therefore 
(see (jnm) also ^^-^ e C"(M3(^'-i)) f^j. q, ^ (q, 1). Clearly, ^ G C°°{R^ \ 

{0}) C C"(M3 \ {0}), and 1^2 ^ (^) = 0, due to the anti-symmetry of the function 
Therefore, all assumptions of Theorem 12.61 are fuUfilled and it follows that 

Secondly, for v^jj^^^i, we make an orthogonal change of coordinates: a = -^{xj — 
Xk),b = -^{xj + Xk), the other coordinates remaining unchanged. Due to the 
specific definition of j, this brings us to a setup exactly as the one above for 

Uj^i. Since the orthogonal change of coordinates does not change the regularity, the 
conclusion follows as before. 

This finishes the proof of Lemma in □ 



Proof of Lemma Vd.dV First, note that the function Gj — j — j satisfies Gj G 
C°(M3W) fQj. a e (0, 1), and 

Gj{xi, . . . ,Xj-i,Xj = 0,0:^+1, . . . ,xn) = 

for all (xi, . . . ,Xj^i,Xj+i, . . . ,xn) G M^^^""^^. 
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Therefore, due to Lemma [2.91 



f • - e C"(M3^) for all a e (0, 1). 
2 \Xj\ 



Secondly, for the function 

1 X.j - Xk 



{(*j-^fe)-<f0,fe)}, 



4 \xj - Xk\ 

the same orthogonal change of coordinates as in the proof of Lemma [3.41 brings us 
in the same situation as the above, again due to the specific definition of *i'(j,fc)- 
The conclusion follows as above. 

This finishes the proof of Lemma □ 

This finishes the proof of Theorem ll.il □ 

3.2. Proof of Theorem 11.41 By Remark 1 1.5 1 is suffices to prove that (|1.20|) holds. 
We proceed similarly to the proof of Theorem ll.il but here we need to estimate 
carefully all the involved quantities uniformly (i.e., independently of xq G R"^^). 
For notational simplicity, we will prove ifL^ only in the case R' = 2R. 
For the proof we need the following regularised version of Lemma 13.11 



Lemma 3.1'. There exists a function Gs.cut : K'^ ^ R, Gs^cut e (^^■^(R'^^), such 
that the function 

Ks^cuti^) {xj-Xk)x{\xj\)x{\xk\)Hx^j+Xk) 

l<j<k<N 

(3.23) + G3,cut(x) = i^3,cut(x) + G3^cut(x) 

(for cut, see (|1.17|) ) solves the equation 

with F2.cut as defined in (|1.16|l and Tcut G G"(R'^"'^) for all a G (0, 1). Furthermore, 
Ga^cut can be chosen such that for all p > the following estimate holds: 

(3-24) l|G3,cut||ci.i(-B3Jv(xo,p)) + lkcut||c°(B3N(xo,p)) ^ C*, 

for some constant C = C{p) > independent of G M^^. 

Proof: The proof of Lemma 3.1' is analogous to that of Lemma [3.11 Instead of 
^,K,v we will use functions /Xcut,«^cut and Vcut to be defined presently. With x 
being the function defined in (|1.14l) we define 

(3.25) iicAx) = x{\A)Kx) ^ x{\x\)\A\ 

(3.26) Kcut(a;, y) = x{\x\)x{\y\)i^{x, y) 

-Jx(3M)(i-x(N))(M^^ 

-Jx(3N)(i-x(M))(N^^ 

2 — TT 

= Xi\x\)xi\y\)^:^{x ■ v) ln(a;^ + y^) + Ki,cut(a;, y). 
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(Note that Ki^cut{x,y) ^ xi\x\)x{\y\)'^i{x7y))- Let I'cut be as in Lemma [B.21 we 
then have 

(3.27) Ai/eut = 73 + h^, 

lkcut||ci-i(B9((a;o,yo,zo),p)) + \\^'^\\c"'(Bai{xo.,vo,zo).p)) < 

with 73 as in H3.5I) and with C independent of (xq, Uo, zq) G and p > 0. 
For /icut, note that 

(3.28) A^cut = AkI' + A(Mcut - 

= 6-A{{l-x{\x\})\x\') = 6-h^, 

where obviously, 

(3-29) ||A*cut||ci.i(B3(:co,p)) + ||^p|lc°(B3(xo,p)) - 

with C independent of xq E M"^ and p > 0. 

For Kcut, using Ak = 72 (see (|IS1) and JISIl), that Aj,(|y|2||^) = 4||i|||, and 
the support properties of X; we have that 

A^ceut = 72 - {1 - X{\x\)x{\y\)}{l - X(3|x|) - x(3|2;|))72 

- {x(3|2/|)(l - X(|:^l)) + X(3|x|)(l - x(|y|)) }(72 + 

+ Ri + R2 + R3, 



X ■ y 



(3.30) = 72 -H72-e(72 + ^n4)+^i+^2 + i?3, 

where 



i?i = x(l2/|)KA,x(|a;|) + x(|y|)2V,x(kl) • V,At 
+ x{\x\)KAyx{\y\) + x{\x\)2Wyx{\y\) ■ ^yn, 

R2^-\xm)\y?r]-^.{a-xiM))^) 

4 \y\ \x\' 

-i(A,x(3M))M^^-((i-x(N)^) 
-^(v,x(3M))-v,(^|yp(i-x(N))), 

and where i?3 is R2 with x and ?/ interchanged. 

Using that k £ C^^"(M^) for all a £ (0, 1), and the support properties of Xi it is 
easily seen that 

(3-31) ll^jllc°=(S6((2;o,yo),p)) < C*' 

with a constant C independent of {xq, j/q)) G and p > 0. 
Since for aU {x,y) £ K.^, 



IV72I < 6\/2f-i- + -i- 



V 72 ' 



^ 8\/2 



\y\J ' |a:||2/|'^ N - 2/| 

we get, using the support properties of Ti. and Q, that 

||HV72||l~(m«)<C , ||c;v(72 + ^)|L^(„.)<C. 
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Again using the support properties of 7i and Q, this imphes that 

(3-32) \\'H'^2\\c<>-HBe{(xo,yo),p)) ^ C", 

with a constant C independent of (a;o, yo) € and p > 0. 
From 12201), (ESU, and we get 

(3.33) AKcut = 72 + , ||^K||c°(S6((a;o,ao):P)) ^ C": 

with a constant C independent of {xo,yo) G and p > 0. Note that (see (|3.26|) 
and ifO)) ') 

'^i,cut(x,y)=x(NI)x(|y|)((2:'+y')G„,(||^)) , G., G ^^'^(SS). 
Therefore, due to the compact support of x, 

(3-34) l|Kl,cut||ci'i(B6((xo,yo),p)) ^ 

with C independent of {xq, yo) E M.^ and p > 0. 
Observe that 

(3.35) = |VF2,eutP + V(F2 - F2,cut) • V(F2 + i^2,cut) 

and that 

V(F2-F2,cut)-V(F2+F2 

.cut ) 

N 

= J2 ^o{F2 - F2,cut) • yj{F2 + i^2,cut) 



N 



j = l ' l<3<k<N ' 



where 



b,=--[l + x{\xj\) + x'{\xM^3\}{(^j+ E 
(j,fe) -|l + x(|a;j - a^fcl) +X'(|a;j - Xk\)\xj -Xfeljx 



Af A 
l=l,l^j l=l,l^k 

= ^{(1-X(I^J -2;fc|)) -x'{\xj -Xk\)\xj -2:fe|||^=^— 

Clearly (using the support properties of x), for all (3 £ N'^^, 
(3.36) ||5^6,|U»(K3«) + ||9'36(,-fe)|Uoe(R3«) < C{f3). 
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Gi, 



2 -^J 



j=i ■ 
i E ■ ( 







72 


\Xj - Xkl-' 



X{\xj ~Xk\). 



l<j<k<N 

Then, due to (|3.36|) and the support properties of x, 
(3-37) l|Gi,cut||ci.i(-B3Jv(xo.p)) < C") 

for some constant C ~ C{p) > independent of xq G M'^^. 
Using A(|.,f ^) = 4^ and A(|^P|^) = 



we see 



(3.38) 
with 



i?, 



, N 3 

j"=i 1=1 ' 

+ i E ^(xCNj -a:fc|)6(j,fc)) • ( 



l<j<*:<Af 
3 



Xj Xk 



V2 



\Xj -Xk\ 



j<k i=l 



Xj Xk 



V2 



■ Xji Xi^^i 
\X, - Xk\ 



N 



E^^-^(i-^(i-^-i)) 

Xj X]^ 



+ E hj. 



\Xj -Xk\ 



From 
(3.39) 



l<j<k<N 

and the support properties of x? we see that 

II^IIco.i(B3n(xo,p)) < G, 



for some constant C independent of Xq G and p > 0. 
Define 



(3.40) 



G3,cut — Gi^cut + G2 .cut 



28 S. FOURNAIS AND M. & T. HOFFMANN-OSTENHOF AND T. 0. S0RENSEN 



G2,cut — Acut + '^l.cut + Vcut, 



1/ ^ Z2 1 
Acut(x) -t( XI -T-^^<^^i'y^]) + XI T^Mcut(a^i - Xk) 



with 



6 4 ^ ^ ^ 16' 

j=l l<j<k<N 

Kl,cut(x) = — X l^l.cutiXj,Xk), 

l<j<k<N 

i>cut(x) = X i^cut(2:j,Xfe,a;;). 

l<j<fc</<Af 

Then, with Ks^cut defined as in we have, using ljT77|l . ijX^ . ijT^ . 

(E3SJ), 1231 

(3.41) - -|VF2,cut|' + rcut, 
where, due to jT^Zjl . (pQ-^ . ^233, (1223, 

(3.42) ||r-cut|ic"(B3„(xo,p)) <C, 

for some constant C — C{p) > independent of xq G M'^^. 
Also, using (12301), 12!ni, (ESni), and (|2!21, 

(3.43) l|G3,cut||ci-i(B3Jv(xo,p)) < C", 

for some constant C independent of xq G M"^^ and p > 0. Now, 13.24|l follows from 
(pn^ and |j2321l- This finishes the proof of Lemma 3.1'. □ 

Let A'a^cut be the function constructed in Lemma 3.1' above. Define (see H3.23|l . 

(trra . and srm ) 

(3.44) C3,cut = e-^=-'-^^'-'V = e-°^'-'<^3,cut. 
Since for all p > (using Lemma 3.1') 

11^3, cut - ^3,cut||ci'i(B3iv(xo,p)) = II Ga^cut || Ci'i (B3]v(xo,p)) 

is bounded independently of xq, to prove (|1.2UI) is equivalent to proving 

(3.45) ||C3,cut||ci'i(B3iv(xo,i?^)) < C'(-R)IIC3,cut||L~(B3iv(xo,2i?,)). 

Using that C3,cut — e^'~^^-'""'4>3,cut, the estimate H3.24|l (twice), and the bound H1.21|l . 
we get, for all < p < p' , 

(3-46) IIC3,cut||ci-°(B3jv(xo,p)) < C'IIC3,cut||L°°(B3jv(xo,p'))' 

with C = C{p,p'). Proving (|3.45|) is improving H3.46|l to a = 1. 
The function C3,cut satisfies the equation 

AC3,cut + 2(VF2,cut + Vi^a^cut) • VC3,cut 

+ (AF2,cut + Aifa,cut + |VF2,cut + Vifa,cut|' + {E - l^))C3,cut = 0. 
We can rewrite this as 

(3.47) ACa.cut + 2VF2,cut • (VCa.cut + C3,cutVX3,cut) 

+ ''l.cut • VC3,cut + '"2,cutC3,cut = 0, 



SHARP REGULARITY FOR WAVE FUNCTIONS 



29 



with (since AF2 = V and Ai^a^cut = -|Vi^2 .cut I ~t~ ?"cut ) 

'"2, cut = ,cut ~t~ ''"cut ~l~ .cut 

= A(F2,cut - F2) + rcut + |Vif3,cutP + 

By the construction of F2 and i^2,cut (see (|1.9ll . (|1.14ll . and (ILIGI) ) it is clear that 
for all p > 

||A(i^2,cut - -F2)||c°(S3jv(xo,p)) < 

with C = C(p) independent of xq G M'^^. Due to Lemma 3.1' (see also H1.14|l V 
Vi^s^cut is C", and we have for all p > 

(3.48) l|Vif3.cut||c=(B3„(xo,p)) <C, 

with C = C{p) independent of xq G K.'^^. This, together with 13.24|l . means that 
(3-49) lkj,cut||c"(S3«(xo,p)) < C*, j = l,2, 

where C — C{p) is independent of xq G M.^^ . 

In order to finish the proof, we introduce a localisation. Let / iM^M, 0</< 
1, be decreasing and such that f{t) = 1 for t < and f{t) = for i > 1, and define, 
for p > 0, A > 1, 

(3.50) 0ix) ^ 0,A'^) ^ /(^(^ - 1)). 

(So 0{x) = 1 on B3n{xo,p) and 9{x) — outside B3n{xq, Xp)). 

Clearly the derivatives of 9 are bounded independently of xq. Below, all constants 
C = C{p) also depend on A > 1; we omit this dependence in the notation. On the 
set B3n{-Xo, p)), ^Cs.cut satisfies the following equation: 

(3.51) A(0C3,cut) + 2Vi^2,cut • (V(0C3,cut) + (0C3.cut)Vif3^cut) 

+ ri,cut ■ V(6IC3,cut) + 7'2,cut(6'C3,cut) = 0. 

Using (|3.51() we will deduce that 

(3.52) l|6'_R,y2 C3,cut||ci.i(S3iv(xo,fl)) < C'(i?)||C3,cut||L°°(B3„(xo,2fl)), 

from which H3.45|l clearly follows (since 9 = 1 on B3n{xq, R)). To prove Theo- 
rem 11.41 it therefore remains to prove H3.52|l . 

Proof of ()3.52|) : Let ^'j,i,cut be defined as ^j^i, was in (|3.12|) but with C3i^3 
replaced by 61(3, cut, ^i^3,cut, that is {j G {1, . . . , N}, i G {1, 2, 3}), 

/Q rn\ ,T, „f^(^'C3,cut) .dKscut 

(3.53) *J,i,cut = 2 — (- 2(6'C3,cut)-^ • 

(Here, 9 = 9^^^). We define *j,i,cut, '^(j.k),i,cnt analogously to ^{j^k)A defined 
in and (|XT^ . Using and we get that for aU < p < p', 

(3.54) ll^iAcut||c.o(i33„(xo,p)) - C'(p)ll^'C3.cut||ci-(i33«(xo,p)) 

< C{p, p', -R)||C3,cut||L°°(B3„(xo,p'))- 

We then have the following result, similar to Lemma [3.41 



30 S. FOURNAIS AND M. & T. HOFFMANN-OSTENHOF AND T. 0. S0RENSEN 



Lemma 3.4'. Let Uj,i.cut, i'(j,fc),i,cut be the solutions to the equations ()3.16(l . (|3.17|l 

(with ^{j,k).i replaced by 5'j,i,cut; ^{j.k).i,cut) given by the Newton potential on 
B3n{xo,V2R). 

Then, for all p < V2R < p' , there exists a constant C = C{p, p' , R) (independent 
o/xo e R^^) such that 



(3.55) 
(3.56) 



Ikj.illci-M-Bsjvlxo.p)) < C'IIC3,cut||L=°(B3„(xo,p'))' 
ll'f(i,fc),illci.i(B3N(xo,p)) < C'IIC3,cut|lL==(B3„(xo,p'))- 



Proof: Using Theorem l2 . 61 and Remark l2 . 71 l(iv|) and (Q, we get the a priori estimate 



||%-^,cut||ci.i(S3iv(xo,p)) < C^SUp 



I* 



sup 



j,i,cut 1 1 c° (ttsjv -3-B3JV (xo ,y2-R.)) 



hi' 



ir3JV-3S3jv(xo,y2-R)) 



j.i, cut 



Xi 



(3.57) 

Using and lpnH|l we have 

||*j,i,cut|lc°(7r3„-3B3Jv(xo,V2fl)) - li*J'».™t||(7°((7r3]v-3B3Jv(xo,\/2i?.))xR3) 

- '^II^C3,cut|lci,»((,r3jv-3-B3iv(xo.V2fl))xR3)- 

This, the compact support of 9, and H3.57|l imphes the estimate 

(3-58) ll"i:i,cut||ci.i(B3Jv(xo,p)) < C" IIC3,cut|lci."(B3iv(xo,V2-R))- 

Combining H3.58|) and H3.46|l . we arrive at H3.55|l . This finishes the proof of the 
estimate (|3.55|) for Uj,i,cut- 

The analogous estimate H3.56II for V(j^k).i,cut is proved in the same manner using 
the same coordinate transformation as in the proof of Lemma l3.4l (see also the proof 
of Lemma 3.5' below). We omit the details. □ 



Lemma 3.5'. Let 5'j,i.cut be defined by (|3.53l) and let ^'j,i,cut o-rid ^(j^k).i,cut be 
defined by 13.14|l and 13.15|l (with i replaced by "^j^i.cut)- Then the functions 
defined by H3.18|l and H3.19|l (again, with an extra index 'cut belong to C"(U.^^) 
for all a G (0, 1). Furthermore, for any p < \/2R < p' , their -norms on the ball 
-B3Af(xo,p) are bounded by 

(3-59) C'IIC3,cut||L°=(S3iv(xo,p')) 

with C — C(p, p' , R) independent of xq e R"^^. 

Proof: That the functions belong to C"(R^^) for all a e (0, 1) follows like in the 
proof of Lemma I3. 51 

To prove the bounds on the norms it suffices, by Lemma 12.91 and the triangle 
inequality, to prove them for 

||*fe,»,cut||c.„(B3iv(xo,p)) ™^ ll*(j''=).^™t)|lc°(i33«(xo,p))- 

For '^k,i,cut, the estimate follows from ()3.54|l . 
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To bound ^(j^k),i,cut, denote by tj^k '■ K K the linear transformation (see 
also (EIEl)), 

(xi , . . . , , 2 ('^j \ -^j+l-j ■ • • 1; 2^*^^ -^/i;)! •^k+li ■ ■ ■ 7 

so that 

Then, since |tj_fc(z)| < |z|, 



|^(j,fc),»,cut(x) - ^0-,fc),»,cut(y)| ^ |^3,»,cut(^j,fc(x)) - ^j,i,cutfe,fc(y))| 



(3.60) 



|x-y|" - |t,,fc(x)-t,-fc(y)|" 

|^'fe,i,cut(ij,fc(x)) - 5'fc,i,cut(ij,fc(y))| 



|i,,fe(x)-ij,fc(y)|" 

Due to the localisation 9 in the definition of ^'^^i^cut (see (|3.53|) '). both of the terms 
on the RHS of (|3.()U|) are bounded by 

cut|lc'i,o,(B3iv(xo,V2-R)) ■ 

The bound H3.59|l for $(j,fc),i,cut now follows using H3.46|l . This finishes the proof of 
the bound H3.59|) for the functions (*j^i,cut - *fc,i,cut) - ^(j,k),i,cut- 

The proof for the functions ^'j,i,cut — *I'j_,;.cut is similar (see also the proof of 
Lemma 3.4' above), so we omit the details. □ 

To finish the proof of Theorem 11.41 define [/cut analogously to H3.20|l . using the 
functions Uj^i^cut, ''^(j,fc),i,cut from Lemma 3.4'. Then, by Lemma 3.4', for any p < 
V2R < p', 

^ Z 1 - 

j = l ' ' l<]<k<N ' ' 

(3.62) l|C^cut||ci-i(-B3iv(xo.p)) C'IIC3,cut||L~(B3Jv(xo,p'))- 

Define [9 ^ 6^^^) 

(3.63) Wcut = 0C3,cut - C/cut, 

then, using (jXKT|) . l|T5S|l . l|TFT|l . and the form of VF2 (see (lOl '). we get the 

following equation for Wcut^ 

Z ^ X 



l<j<fc<7V ' ' 



AT 



+ ^ (F2 - i^2,cut) -^j.cut 

i=i 

(3.64) ~ |ri,cut • V(6'C3,cut) + r2,cut(6'C3,cut)} = A. 

Here, A belongs to C" for aU a e (0, 1), and, for all p < \/2R < p' , 

(3-65) ||A||c°(B3n(xo,p)) < C'IIC3,cut||L=°(B3N(xo,p')) 
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with C = C{p, p' , H) independent of xq G M.^^ . For the first two terms in H3.64|l this 
follows from Lemma 3.5'. For the third term, it follows using the form of F2 — F2,cut 
(see (fT^ . (11331, and ifTTCIl l and For the last term we use and 

By Proposition 12.41 this means that Wcut belongs to C^'", and we have the 
estimate 

(3.66) |lWcut||ci>i(B3iv(xo,fl)) < llW^cut|lc2,°(S3„(xo,fl)) 

< C{R){\\Wcnt\\Lo.(^BsN(^o,^R)) + ll^llc°(B3iv(xo.^fl)))- 

Using H3.63|l , the triangle inequality, and then (|3.62|) (with p — \/2R and p' — 2R) , 
we have 

||Wcut|li^(B3„(xo, 3/2fl)) - C'(i?)||C3,cut||L~(S3jv(xo,2i?))- 

This, (EiSEl), and IpTHKj) with p = ^/2R and p' = 2R, gives the estimate 

(3.67) l|Wcut||ci'i(B3]v(xo,-R)) < C'(^)IIC3,cut||L°°(S3jv(xo.2i?.))- 

Using 6'C3,cut = Wcut + C/cut, (with p ^ R and p' = 2i?) and 1(233, the 

estimate (|3.52f) follows. □ 

This finishes the proof of Theorem 11.41 □ 

Appendix A. Construction of the function k 

In this appendix we construct the function k, that gives rise to the terms of 
order ln(r) in the function solving = — |Vi^2p (see the previous section, 
Remark 13.31 in particular). Therefore, k is responsable for the '"-singularities in 
the wavefunction ip. 

More precisely, we prove the following: 

Lemma A.l. Let the function 72 : — > K he given by 



(A.l) '^2(2;,y) = (-r-r--r-r •] r , x,y £ 

\\x\ \y\J \x-y\ 

: exists a function k :MP ^ R of the form 
«(x, y)^^^ix-y) log{x' + y^) + [x' + y^) (i^) 



Then there exists a function k : M'' ^ M of the form 

^— ^ {x ■ y) \og{x^ + y^) + [x' + y^) 
37r ^\{x,y)\ 

(A.2) =k{x,y) + Ki{x,y) , C £ C^^\^'') 

satisfying Ak — 72 . 

Remark A.2. Note that by Lemma WHA ki e C^'^I*^)- 

Proof RecaU that t)f^ = Ran(7'f ^) is given by the linear span of the harmonic, 
homogeneous polynomials of degree 2 in restricted to 
By Lemma FC. II we have that 

/T,(6) ^/ N 16(2 -7r) 

(^2 72)(-) = ci-^ , cr = ^^, 

where = + y'^,u! = {x, y)/r e Let k{x,y) = j^{x ■ y) \og{x'^ + y^). Then 
(A, + Ay)k{x,y) = ci^-^ = (T'f 72)(ra>). 
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Letting Ki = K—k this reduces the problem (of finding k such that (Ax + Aj,)k = 72) 
to finding ki such that 

(A.3) (A^ + Ay)K,i = 72 

with 

(^•4) 72 = 72 - ci . 



Due to the above, 



(p(^H2)(rc.) = 0. 

Therefore, by Proposition l2.8l there exists a solution ki to ljA.3|) such that Ki(r(jj) = 
r^GM. with e Ci>"(§5) for all a e (0, 1). 

To verify ifOjl we need to prove that in fact G^i G C^'^(§^). We will do this by 
proving that m e C^^^{M.^ \ {0}), since then G^, = Ki/r^ e C^^^S'^). 

To prove ki e C1'^(R6\{0}), we analyze the equation lfO|) for Ki in the vicinity 
of singular points of the function 72 on the sphere There are two types of singular 
points: (a) {xo,a;o) G S^, (b) {0,yo) € (resp. {xo,0) G §^). The function ki is 
C°° in a neighbourhood of all other points on §^ due to Proposition l2.4l fsince. for 
r > 0, 72 is C°° away from points of type (a) and (b), see l)A.l(l and l|A.4(l 'l. 

(a) : Let Ua C be a neighbourhood of a point {xq,xo) £ (i.e. 2\xo\^ = 1) 
such that for some c > 0, > c, \y\ > c for {x,y) G Ua- Choose new coordinates: 
Let 

{xi,X2) =t{x,y) = (x-y,x + y). 

Then 

/ ,_iw ^ Xi f Xi- X2 Xi+ X2 \ _ Xi 

with Ga e C°°(t(C/a)). Since GaiO,X2) = for X2 7^ (that is, for a; = y 7^ 
in the original coordinates), we have, by Lemma f2. 91 that 72 o t^^ g C*^'^ (t{Ua)) , 
and therefore 72 G G°'\Ua) C C"(J7a) for aU a G (0, 1). Since (a; • y)/{x^ + y^) G 
C°°{Ua), we have (see (EIJ) 72 G C"(C/a) for all a G (0, 1). By Proposition ITl we 
get from (fO)l that ki G C^-" ([/„). 

(b) : Let Ut C be a neighbourhood of a point (0,?;o) G §^ (i.e. \yo\ ~ 1) such 
that for some c > 0, \y\ > c,\x — y\ > c for (a;, y) G Ut- Then 

y \ a; - y 



X y X f y y — X \ y 



\x\ \y\ \x\ \\y\ \y-x\y \y 



Note that 



and that 



G C^iUb) 



\x\ Myl \y-x\) ~ \x\ ^ ' ' 
with Gb G C°°{Ub), G6(0, y) = for 2/ 7^ 0. Therefore, by Lemma 1^ and (TOt . 

72(x,y) - ( - ^ • ^) G GO'i(C/b) C G"([/b) for all a G (0, 1). 
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Let K2 be such that 



\x\ \y\ 

The existence of such a function is ensured by Theorem 12. 61 since y for {x,y) S 
C/b, and V'^^ (]§[)" ^ '^^^ anti-symmetry of 

Then (see (|A.3(l 'l kz — k,i~ K2 solves 

(A, + A^)«3 = 72(a:, y) - ( - ^ • ^) e C"(C/b) for aU a e (0, 1), 

so by ehiptic regularity K3 € C^'"(?7b) C C^'^(?7b). Since K2 G C^'^(C/b), this proves 
Ki = K2 + G C^'^ ([/(,). Together with ki G C^'"(?7a) from above, this implies 
G^, = Ki/r2 e Ci'i(§5), and so ki = G Ci^i(RS). 

This finishes the proof of the existence of k solving H3.6() . and having the form 
(fO^ . with G = G.,. □ 

Appendix B. Construction of the function v 
In this appendix wc construct a function solving H3.7|l . 

Lemma B.l. There exists a solution v = v{x,y, z) to the equation 1)3.7(1 satisfying 

(i) v is invariant under cyclic permutation, i.e., v(x,y,z) = (i/ o cr)(x,y, z) for 
all x,y, z £ 'W? , where a{x, y, z) = (z, x, y). 

(ii) z/G Ci'i(K^). 

The idea is to change coordinates, to the centre-of-mass frame for {x,y,z). In 
these new coordinates, the problem of solving ((3.7|l turns out to reduce to a problem 
in 6 variables only. By an extra symmetry of the function 73 (see (|3.5|l V namely 
permutation of the three electron-coordinates x,y, and z, the logarithmic term 
that occured in the function k (see (|A.2ll 'l does not occur here. This is because the 
projection on t)^^ of 73 (the function that 73 transforms into in the new coordinates, 
see (IB.2I1 below) vanishes, due to this extra symmetry. 

Proof: Make the following change of coordinates (each entry below is a diagonal 
3 X 3-matrix with the listed number in the diagonal; we will use this notation 
repeatedly; here, x,y, z G M.^) 





1 





2 








1 










V3 


V2 


V6 









X2 




\ -3 / 



(B.l) \ y =r 



Then 

(B.2) (73or)(xi,X2,a;3) = 

X2 X2 + V^X3 X2 X2 — V^X3 X2 + V^Xs X2 — V^X3 



\x2\ \X2 + V3x3\ \x2\ \x2~V3x3\ \x2 + V^xal |a;2-V3x3| 
= 73(a;i,2:2,a:3). 

That 73 is independent of xi is the fact that 73 only depends on the inter-electron 
coordinates {x~y, y—z, z—x respectively), and not on the centre-of-mass coordinate 
{xcM ^ ^{x + y + z) ^ Xl). 
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The function 73 is invariant under cyclic permutation of the electron-coordinates 
X, y and z, that is, (73 o a) {x, y, z) — 73(0;, y, z) for all x,y, z £ M"^ with a{x, y, z) — 
{z,x,y). This gives that 

(B.3) (73 o n){xi,X2,X3) = 73(2:1, a;2,X3) for aU xi,X2,X3 S R^, 

with TZ the orthogonal transformation given by 7?, = o a oT , that is by the 
9 X 9-matrix (again, each entry is a diagonal 3 x 3-matrix) 

/I 

7^= cos(f^) sin(^) 
V "Sin(f ) cos(f ) 

Note that 7?. is a rotation of (3:2, x^) by ^ around xi (all in R^), that is, 7i? ~ Ig, 
where /g is the identity on R^. 
Define the function 73 by 

(B.4) 73(2:2,2:3) = 73(2:1, 2:2,2:3) , (0:2,2:3) G R*^ 

(since 73 is independent of xi, this is well defined). Then, due to (jB.3p . 

(B.5) (73 o^) (2:2, 2:3) = 73(2:2,2:3) for aU a;2,a::3 G R^, 

with (each entry still being a diagonal 3 x 3-matrix) 




cos(2f) sin(^) 
-sin(f) cos(f) 



(B.6) n = 

3"^ ^"=V — 

Observe that if D — D{x2,X3) solves (for 73, see HB.2|) and (|B.4ll 1 
(B.7) (A,, + A,3)z> = 73, 

then trivially the function D defined by i>(xi, X2, X3) = D{x2,X3) solves 

Since T is orthogonal, the function v = v o T^^ will then solve (recall that 73 = 
73 oT) (^Ax + Ay + AzVu = 73, that is, H3.7|) . The problem of solving H3.7|) therefore 
reduces to solving (|B.7p . 

Observe next that (see HB.2II and (|B.4|) ') 

73(0x2,02:3) = 73(2:2,2:3) for aU O e 5*0(3), 0:2, 2:3 G R^ 

This and ||R5|| gives, by ^ of Lemma El that vf^^^ = . Therefore, by 
Proposition 12. 81 there exists a solution v to (|B.7|I with 

-/ N ( 2 , 2-,r< ( (2^2,2:3) 
1^(2:2,2:3) = (2:2 + 0:3)65 



Vl(2:2, 2:3)1, 

gC1'"(S5) for all a G (0,1). 

We proceed to prove that in fact Gp G G^^"(S^) for all a G (0, 1). We do this by 
showing that v G G^'"(R^ \ {0}), using IIB.7|) and elliptic regularity (Proposition 

El. 

Note that there are two kinds of singular points of 73 on §^ : (a) 2:2 (and so 
2:3 7^ 0), (b) X2 = V^X3 (and so 2:2 7^ 7^ 2:3) (resp. 2:2 — —^32:3). The function 
V (and therefore, Gp) is G°° in a neighbourhood of all other points on S'"^ due to 
elliptic regularity (Proposition 12. 4|1 . 
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(a): Let Ua C be a neighbourhood of a point (0, x^) E §^ (i.e., ^ 0), such 
that for some c > 0, \x2 + VSx^] > c, \x2 — V^x^] > c for {x2,X3) £ Ua- Note that 

_ X2 I X2 + VSxs X2 - VSX3 \ 

13(X2,X3) = -. r 



(B.8) 
Write 



X2\ \\X2 + V3X3\ \X2-V3X3\ 

X2 + \/3a;3 X2 — Vixz 

\x2 + \/iX3\ \x2-\/iX3\ 
X2 I X2 + V^X3 , X2 - V3X3 \ X2 



j- • Gaix2,X3) 



\X2\ \\X2 + V^X3\ \X2~V3X3\J \X2 

where Ga £ C°°{Ua), GaiO,X3) = 0. Furthermore, 

X2 + V3X3 X2 - V3X3 ^ (jooijj -J 
\X2 + Vix3\ \x2~ViX3\ 

Therefore, due to Lemma El 73 £ C°'i(f/a) C G"{Ua) for all a £ (0, 1), and so, 
by (|B.7p and elhptic regularity (Proposition EH), v £ C'^^°'{Ua). 

(b): Let Ub be a neighbourhood of a point (0:2, 2^3) £ with X2 ~ x/Sxjj (i.e., 
2^2 7^ 7^ ^3); such that for some c > 0, |a;2| > c, \x2 + V^X3\ > c for {x2, X3) £ Ub- 
Choose new coordinates: Let 

{u,v) = t{x2,X3) = {X2 - y/3x3,X2 + Vixs). 

Then 

U / U + V V \ U + V V 



(- N u f U + V V \ 

(73 OT )(u,w = — • ■ r--r-r 

\u\ \\u + v\ \v\J 

We proceed as above. Write 



+ |w| 



u\ \ m + v\ \v\ / \u 



u+v V \ u s 
Gb(u,v) 



where Gb £ C°°(t(C/6)) (since v ^ 0,u + v ^ in T{Ub)), Gb{0,v) = for d ^ 0. 
Furthermore, 

" + " ^ £C-[Ub). 



\u + v\ \v\ 

Lemma El implies that 73 o £ G°^^{T{Ub)), and so 73 € G^^^{Ub) C C"(J7fc) 
for all a £ (0,1). By (|B.7|) and elliptic regularity (Proposition I2.4|l follows that 
u£C^''^{Ub). 

Singular points of the form x^ — — \/3a^3 are treated analogously. 
From the above follows that v £ C2'"(M6 \ {q}), and therefore Gp £ C^^°'{S'^), 
for all a £ (0,1). 

This finishes the construction of a function D £ C^'^(M®) that solves (|B.7|) . and 
has the form 

(B.9) p(a;2,X3) = (x2 + x2)G,(^||^), 

Gp eG2'"(S5) for aU a £ (0,1). 

As discussed above 17 defines a function 1/ solving the equation H3.7|l . Clearly, 
since V £ G^'^(]R^), we get v £ C^-^{R^). The solution 1/ constructed in this 
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manner does not necessarily satisfy the invariance property (jlj. In order to force 
this invariance, we consider 



1 ^ 



Since the Laplace operator commutes with a, and 73 is invariant under cr, t'sym 
satisfies the conclusion of Lemma IB. II □ 



With the notation from the proof of Lemma [B. II we define 

i7cut(a;2,a;3) = xi^l + xl)V{x2,X3), 

with X in (|1.14|l . and i'cut(a;i, X2, 0:3) = X'cut(a;2, 2:3) (as already defined). As 
discussed above (for v) the function Vcut defines a function Vcut = i^cut o : 
^ R (by the linear transformation T in (|B.l|l '). We then get: 



Lemma B.2. The function Vcut satisfies 

with 73 as in (|3.5(l and h G C"(R^) for all a e (0,1). Furthermore, we have the 
estimate 

(B.IO) \\l^cut\\c^-^Bs{{xo,yo,zo),R)) + II ^11 C°(S9((xo,yo,^o).fl)) ^ C*' 

with C independent of {xq, yo, zq) G and R > 0. 



Proof: We calculate, using ljB.7(l . 

(A^^ + A^^ + Aa;3)i>cut = {Ax2 + Ax.^)Vcut = AFcut 

= 73 + { [Axp + 2Vx • V77} - (1 - x)% 

= 73 + /l. 

Using l|B.8|l and (|B.9|I we see that the term in {•} is C" and has compact support. 
The function (1 — x)l3 is (this was proved in the proof of Lemma [B.1|I and 
homogeneous of degree zero outside -66(0, 2). Therefore, 

\\Hc''(Ba({x'>,x'lx°),R)) < C", 

with C independent of {xi,X2,x'^) G R^ and R> 0. Since x has compact support, 
andFe Ci'Hl^^)> we have 

\\'^cut\\c^''^iBa({x°,x°,x°),R)) < d, 

with C independent of [x^, X2, a;^ ) e R9 and i? > 0. 

Since T is an orthogonal transformation, IjB.lOjl follows. This finishes the proof 
of the lemma. □ 
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Appendix C. Computation of 7^2^^72 

In this appendix we compute 7^2^'' 72, the singular part of the two-particle terms 
in IVF2P, see HH.4|I and (|3.5(l . This is Lemma \C1\ below. It follows from general 
results on 7^2^^ rj when has certain symmetry-properties fLemma lC.2(l . The latter 
is also responsable for the non-occurence of terms of order ln(r) (of regularity 
C^'" only) in the function i/ constructed in the previous appendix; see Lemma FB. II 

Lemma C.l. Let 

(C.l) ^,(:,^y) ^ - y) . {x,y)eR'xR^ 

\\x\ \y\J \x-y\ 

Then 

oTT X ~v y 

Proof: This will follow from Lemma FC. 21 and Lemma FC.3l below. Namely, by (P) and 
(O in Lemma IC. 21 we get that, due to symmetry, 

(^2^^72)(2:, y) = ci ^ ^ ^ for some Ci € R, 
X -r y 

that is, only the function x ■ y (restricted to contributes to the projection onto 

\)"2^ of the function 72 in (|C.1|) . That ci = "^^^.^ ^"^ is the result of Lemma |C.3I 

(which is merely two computations). □ 

Lemma C.2. Assume rj E L^(S^) satisfies 

(C.2) r,iOx.,Oy)^fjix,y) 

for all O e 50(3) and almost all {x,y) e §^ C x M^. Let Qi be the orthogonal 
projection ( in L^ ) ) onto 

Span{Pi|g5 , -P2IS5} , 
and Q2 the orthogonal projection onto 

Span{Pi|g.}, 

where Pi{x, y) — x ■ y, P2{x, y) — x"^ ~ y^ , (x, y) G x M^. 
Then 

(i) V'i\ - Qiry. 

(ii) Let r\ satisfy 



(C.3) fl{x,y) = viViX) for almost all {x,y) G S"^ C x 

Then V^^rj — Q2r]. 
(iii) Let TZ be as in ljB.6p . Assume rj satisfies 



(C.4) v{T^{x, y)) = ri{x, y) for almost all {x, y) e §^ C M"* x 
Then V^2\ = 0. 
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Proof of Lemma \C.S\ : Suppose Q is proven then the proofs of Jul) and (Imjl are 
simple: 

Proof of ^ : Due to we only need to prove that 

l{x,y){x^ -y^)duj = 0. 



This follows using the symmetry ljC.3|l of rj (which preserves the measure du of S'^): 

ri{x,y)P{x,y)duj = I- [ f]{x,y){P{y,x) + P{x,y)) du), 

and when P{x,y) — P2(x,y) = x^ — y^, then P{y,x) + P{x,y) — 0. This proves 
©. □ 



Proof of (jnl|l : Using ^ and ljC.4|) it is enough to show that 

P{x,y) + P{1l{x,y)) + P{n\x,y)) = 

when P{x,y) — x ■ y or x^ — y^ (since TZ preserves the measure of S'^). This 
follows by direct calculation. □ 

It remains to prove (0): 

Proof of (0) : Recall that = Ran(-pf ^). Define f)2,m« by 

^2,^nv = Span{/ e I f{Ox,Oy) = f{x,y) for aU O G S0{3)} . 

Note that 7^2^^ ry G i)2,inv because of HC.2|) . We need to prove that 

i}2Mv = Span{Pi|g5 , ^'2ls5}. 

Since every function in l)2,inv can be written as a finite sum of spherical harmonics of 
degree 2 it suffices to consider a real, harmonic polynomial P which is homogeneous 
of degree 2, and which is invariant under the action of S0{3): 

(C.5) P{Ox, Oy) = P{x, y) for aU O E S0{3). 

Identifying P with a quadratic form on M^, there exist real symmetric matrices 
A, B, and C, such that 

(C.6) P{x,y) ^ X ■ Ax + y ■ By + X ■ Cy. 

The condition of harmonicity of P becomes Tr[yl + B] ^ 0. We prove that A, B, 
and C have to be multiples of the identity matrix I3 on M.^. To do so, let us first 
restrict to x = 0. Using (|C.5|I and (|C.6p we get 

yBy = P(0, y) - P(O0, Oy) = Oy ■ BOy, 

for all O G 50(3). Let A be a (real) eigenvalue of B, with corresponding eigenvector 
v. Bv = Xv. Let y be any vector in M'^. Then there exists an Oy G S0{3) such that 
Oyy ~ fiyV for some ^y £ M, and therefore y ■ By = Oyy ■ BOyy = A|jy||^. Since 
this is true for all y G M.^, we get B = A/3. A similar argument (with y — 0, and 
letting X vary) shows that also ^ is a multiple of the identity. Finally, the condition 
of harmonicity, Tt[A + B] = Q, implies that A = —B = —A/3. 
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Finally the term x ■ Cy. This will be treated similarly. Due to the above (see 
(|C.6|I ). X - Cy = P{x,y) - X{y'^ - x'^). Therefore, l|C.5(l implies 

x-Cy = Ox- COy for all O € 50(3). 

By arguments similar to the above, we find that C is also a multiple of the identity 
73. Since P{x,y) — X{x^ ~ y'^) + x ■ Cy, this finishes the proof of iP). □ 

This finishes the proof of Lemma FC. 21 □ 

Lemma C.3. Let Q2 be the orthogonal projection (in L^(§^)j onto 
Span{Pi|g5} , Pi{x,y)=x-y , {x,y) E R"" x , 

and let 
Then 

Proof: Note that, with 



-Pils= IL2(S5) ' Va:2 +2/2' 

we have |li^|lL2(ss) — 1, and so 



(C.8) Q272(w) = r(t^) / Y{u;)-f2{uj)du; 

PI ^2 / ^ils=H72(^)d^i --^T^. 

IL^(s=) j ^ +y 

We need to compute the two integrals in the brackets. 

Since Pi is homogeneous of order 2 and 72 of order (as functions on M^), we 
have 

Pi{x,y)i2{x,y)dxdy = — I Pi|g5(w) 72(w) dw. 
Therefore, 

(C.9) / Piy{cj)^2{io)duj = 8 [ Pi{x,y)^2{x,y)dxdy. 

Js^ J Be (0,1) 

Choose coordinates (|a;|, \y\, \x — y\,^l) for (with Q three necessary angles). Note 
that 

Pi(a;,y)=x.2/ = i(|a;p + |y|2-|x-y|2) , {x,y)e^^-^^ 



X 



and 
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Then (see Hylleraas [H (45d)]; let s = \x\,t = \y\,r ^ \x - y\) 
I Pi{x,yh2{x,y)dxdy ^ dil 



(04) 

1 /-Vl-s^ /-s+t 

{s'^ + f - r'^){s + t){2st~- {s'^ +f - r^)) drdtds 

Jo As-A 



Using (jC.9|l and (jC.lOp this means that 

(C.ll) £ p,|^,(^)^2(^)d^^ ^ j dQ. 

Next, observe that, again due to homogeneity, we have 

?io 



/ {x ■ yfdxdy = — ||Pi 



and so 

(C.12) ll^l|sB||'..s5^ = 10 / [x-yfdxdy 

Since x ■ y = |(|a;p + |yp — 1^; — y\^) we get (using coordinates as above) 
(a; • y)^ dx dy 



\/l-s2 ps+t 



+ - r'^Y srt dr dt ds 



\s~t\ 



dn. 



1280 

This means (see (lUn^ ) that 

II I l|2 ^ f 

||-Pl|s5 ||l2(S5) ^ J 

Now (EH follows from fCl8)l . fCnTjl . and (jHH- This finishes the proof of 
Lemma IC.3I □ 
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